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Motivated by a recent experiment in a tunable graphene analog [L. Tarruell et al, Nature 483, 
302 (2012)], we consider a generalization of the Landau-Zener problem to the case of a quadratic 
crossing between two bands in the vicinity of the merging transition of Dirac cones. The latter 
is described by the so-called universal hamiltonian. In this framework, the inter-band tunneling 
problem depends on two dimensionless parameters: one measures the proximity to the merging 
transition and the other the adiabaticity of the motion. Under the influence of a constant force, the 
probability for a particle to tunnel from the lower to the upper band is computed numerically in the 
whole range of these two parameters and analytically in different limits using (i) the Stiickelberg 
theory for two successive linear band crossings, (ii) diabatic perturbation theory, (iii) adiabatic 
perturbation theory and (iv) a modified Stiickelberg formula. We obtain a complete phase diagram 
and explain the presence of probability oscillations in terms of interferences between two poles in 
the complex time plane. We also compare our results to the above mentioned experiment. 
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I. INTRODUCTION 

The recent surge of interest in bandstructure of un- 
usual crystals spurred from various developments in 
studying condensed matter systems. To name two ex- 
amples, relevant to the present work, are the successful 
isolation of single layer graphene pQ and the fabrication of 
3D topological insulators [3J [3] . In contrast to ordinary 
two-dimensional (2D) crystals, the low-energy hamilto- 
nian involves at least two coupled bands. It gives rise to 
band crossings which, depending on the material param- 
eters, can be gapped or not. The extra degree of freedom 
in the internal space (generically called the pseudospin 
space) offers opportunity for the investigation of new sys- 
tem properties. 

Most recently, hamiltonian engineering with artificial 
crystals - which is a main theme of the active field of 
quantum simulation, see e.g. [4] - provides a complemen- 
tary route to realize coupled-band systems [5HZ]. While 
the low-energy hamiltonian mimics closely that of its 
solid-state counterpart, it is no longer limited to param- 
eters of the actual material. For example, in the cold 
atom experiment performed at ETH Zurich [5J, a Dirac 
cone pair in the bandstructure is brought to merge as 
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a function of laser parameters, thus realizing a Lifshitz 
transition which has never been reached in graphene [5J- 
lllj . The merging is a topological transition in which two 
Dirac cones of opposite Berry phase approach and anni- 
hilate before a gap opens. 

In the ETH Zurich experiment, Bloch oscillations of 
non-interacting fermionic atoms in a honeycomb-like op- 
tical lattice are executed to study the merging transition 
of Dirac cones. As the atom traces out a closed trajec- 
tory in the momentum space, it may tunnel to the second 
band when it comes close to a linear avoided band cross- 
ing (i.e. a Dirac cone), a process known as Landau-Zener 
tunneling. By measuring the transfer probability after 
performing a Bloch cycle, information about the band- 
structure can be extracted with momentum resolution 

In Ref. [13], we presented a tight-binding model that 
reproduces the optical band structure in the parameter 
space of the experiment. Using a low-energy description 
of the tight-binding model known as the universal hamil- 
tonian [11] , we quantitatively reproduce the experimen- 
tal results of Ref. [HJ . In the framework of the universal 
hamiltonian, the inter-band tunneling problem depends 
only on two relevant parameters: the merging gap d - 
which controls the proximity to the transition - and the 
momentum perpendicular to the direction of motion k - 
which controls the adiabaticity of the motion, or in other 
words how far in momentum space the atom is from hit- 
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ting exactly the tip of the Dirac cone. In particular, in 
Ref. [13] we explain the situation where the two Dirac 
cones are hit in succession during a single Bloch oscilla- 
tion (see Fig. [I]) by using a simple approximation, known 
as the Stuckelberg theory EH], m which tunneling 
events are assumed to be independent. The validity of 
this approach is restricted to the gapless phase (d < 0) 
and not too close to the merging transition {d <C —1,—k), 
i.e. the two Dirac cones are well separated. 

The present paper is an extension of our letter [13] and 
focuses on the tunneling problem where the atom encoun- 
ters two Dirac cones in succession. Here, we go beyond 
the independent cone approximation and present a com- 
plete picture of the inter-band transition probability as 
a function of the two parameters d and k. In particular, 
we now access the whole phase diagram, including the 
gapped phase (d > 0) and the transition point (d = 0). 
Our paper is organized as follows. In section II, we formu- 
late the inter-band tunneling problem for the universal 
hamiltonian. In section III, we recall the approximate so- 
lution used in [13] based on the Stuckelberg theory. We 
then present three other analytical approaches: diabatic 
perturbation theory in section IV, adiabatic perturbation 
theory in section V and a modified Stuckelberg formula 
in section VI. In section VII, we present numerical solu- 
tions in the whole parameter space and compare the re- 
sults of the different approaches. Finally in section VIII, 
we compare our results to the ETH Zurich experiment 
before concluding in section IX. 




FIG. 1: (Color online) Energy spectrum in the gapless phase 
(A* < 0): energy E = ±^/(pl/(2m*) + A*) 2 + c%y% as a 
function of momentum p x ~ t and p y ~ k. The distance 
between the two Dirac cones is controlled by the merging gap 
A, oc d. The perpendicular gap c y p y oc k controls how far 
the particle is from hitting the Dirac cones, which are located 
at (t = ±^/|d[, k = 0), directly at their tip. The black lines 
are lines of constant k. 



INTER-BAND TUNNELING IN THE 
UNIVERSAL HAMILTONIAN 



II. 



In the Landau-Zener (LZ) problem [16 - 18 , an avoided 
linear crossing between two bands is considered and the 



probability for a particle (which will be called an elec- 
tron in the following) to tunnel from the lower to the 
upper band under a constant applied force is calculated. 
Landau solved the problem approximately using pertur- 
bation theory and the semiclassical approximation [16] . 
while Zener was able to find the exact solution [17] . For a 
concise modern presentation, see Ref. 18 . Here we con- 
sider such a tunneling problem for the case of a quadratic 
band crossing. The latter occurs close to the merging 
transition of Dirac points [STfTT] and was recently ob- 
served in a cold atom realization of a graphene analog 
[5J. In the gapless phase, the quadratic band crossing 
can be approximated as two successive linear crossings 
(or Dirac cones), which is at the heart of the Stuckelberg 
approach (see below). 

We start from the universal hamiltonian describing the 
vicinity of the merging transition |llj : 
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It depends on three real parameters: an effective mass 
m* > giving the spectrum curvature in the x direc- 
tion, an effective velocity c y > for the y direction 
and a merging gap A* , which is a real number con- 
trolling the distance to the transition EH]. The state 
space is that of an electron moving in a two-dimensional 
plane and carrying a pseudo-spin 1/2 described by the 
Pauli matrices a x , a y , a z . The corresponding spectrum 



is E = ±y (p x /2m* + A*) 2 + c y p y and is plotted in Fig. 

[1] when A* < 0. If the merging gap is negative, the 
spectrum is gapless and contains two Dirac cones at 
(p x = ±y / 2ro7|A7f,Pj, = 0). If it is zero (the merging 
point), the two Dirac cones are on top of each other and 
the spectrum is linear in one direction and quadratic in 



the perpendicular direction E = ±^J (p x /2m«) 2 + c 2 p 2 

[3]. If it is positive, there are no band touching points 
anymore but a true gap 2A* between the two bands. 

We add a constant electric field £ in the x direction 
such that during its motion an electron encounters the 
two Dirac cones in succession [5U] , see Fig. [I] The gauge 
is such that the vector potential A x = —£t and A y = 0. 
Therefore 



H u {t) = 



2m* 



■A, 



(2) 



The force F — e£ is taken to be positive and — e < 
is the electron charge. The hamiltonian commutes with 
p x and p y and therefore the non-trivial dynamics only 
occurs in the internal space of the pseudo-spin 1/2 and 
p x and p y can be taken as c-numbers (conserved quan- 
tities). Shifting the origin of time Ft — p x — > Ft, it is 
now possible to get rid of p x . This hamiltonian defines 
a characteristic energy scale E c h ar = (hF) 2 ^ 3 /(2m*) 1 / 3 , 
and therefore a timescale t c har = h/E c h a r and a length 
scale L c h ar = E^^jF . Energies, times and lengths are 
therefore given in units of these characteristic scales. We 
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then define the dimensionless quantities d = A*/E c h a r 
and k = c y p y / E c h ar and the dimensionless hamiltonian 
H u (t) = [t 2 + d]a x + ka z . 

Performing a unitary transformation in pseudo-spin 
space allows one to rewrite the 2x2 hamiltonian in a 
familiar LZ form. Let (<J x ,a y ,a z ) — ¥ (c z , <r Xl a y ) which 
is realized by the unitary operator U — exp(i^cr ■ n) = 
^(I+ia x +i(Ty+ia z ) where n = (1, 1, l)/y/3. Then H u (t) 
becomes 

H V=( E £i + + (3) 

where Ei(t) = —E2(t) = t 2 + d is a quadratic function 
of time (in contrast to the original LZ problem in which 
Exit) = -E 2 (t) oc t) and H 21 = H 12 = k = const. The 
quantities d and k are the only two relevant dimensionless 
parameters. The first parameter, d, controls the distance 
to the merging transition, which occurs at d — 0. When 
d < there are two Dirac cones (gapless phase) and 
when d > there are no Dirac cones (gapped phase). 
The other parameter, k, controls how far the electron is 
from hitting the Dirac cones directly [H] , see Fig. [TJ it is 
also a measure of the adiabaticity. We call d the merging 
gap and k the perpendicular gap. The orthonormal basis 
{ 1 1) , 1 2) } in which the hamiltonian is written is called the 
diabatic basis. The diabatic spectrum corresponds to a 
negligible k and is simply E\(t) = t 2 + d and E 2 (t) = 
—Ei(t). It is plotted in Fig. |2^a) for negative d. It 
features two band crossings in real time. 
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FIG. 2: (Color online) Energy E as a function of time t when 
d — —2 < 0. (a) Diabats Ei,2 = ±[i 2 + d] are parabolas 
that intercept in real time (continuous blue line is for + sign 
and dashed red line is for — sign), (b) Adiabats E+ t - = 
±^J\t 2 + d] 2 + k 2 (with k = 0.5) do not intercept in real time 
(continuous blue line is for + sign and dashed red line is for 
— sign) but do intercept in complex time. 



The state of the electron at a given time t is described 
by the bispinor \ip(t)). Its time evolution is given by the 
Schrodinger equation i4z\ip) = H(t)\ip). Let us assume 
that initially (t — > -co), the electron is in the lower band 
\ipi~ oo)) ~ |2). Our aim is to compute the probability 
P(k,d) = \(l\ip(oo))\ 2 that it ends in the upper band 
(t — > oo) as a function of the two parameters k and d. 
As the result does not depend on the sign of k, we will 
assume that k > 0, without loss of generality. 

In the following, we mathematically formulate this 
problem in two different bases, namely the diabatic and 
the adiabatic bases. 



A. Diabatic basis 

In the diabatic basis, we write the state at an arbitrary 
time as a function of two complex numbers Ai(t) and 
A 2 (t): 

\if>it)) = A 1 {t)e-^ tdt ' E ^ t ^\l)+A 2 {t)e- l ^ dt ' E2{t ^\2) 

(4) 

The time evolution is governed by the Schrodinger equa- 
tion, which reads: 

A x = -iH l2 A 2 it)e i S tdt ' E ^ t ' ) 

A 2 = -iHl 2 A l it)e- i S tdt ' E ^ t ' ) (5) 

where E\ 2 (t) = E\{t) — E 2 (t). One therefore needs to 
solve this system of coupled equations with the initial 
conditions Ai (— oo) = and A 2 {— oo) = 1 (up to a global 
phase factor). As |Ai(i)| 2 + |A 2 (t)| 2 = 1 at any time, we 
are only interested in finding P = |yli(oo)| 2 . This system 
of two coupled first order differential equations can also 
be written as a single second order differential equation 
for A\ (or A 2 ) alone [17] . 

If the force is large, the motion of the electron is fast 
and k = CyPy^m*) 1 / 3 /(hF) 2/3 < 1 is negligible. This 
is the diabatic or sudden limit. In such a limit, the elec- 
tron stays in the lower state |2) and P — > 0. Indeed, 
when k = 0, the tunneling probability is zero for all d as 
a result of the conservation of the pseudo-spin a z , which 
commutes with the hamiltonian H{t). This may seem 
surprising as it means that even when the two bands 
overlap the probability of interband tunneling is zero. In 
particular, when d — with a single quadratic crossing 
point, the electron does not tunnel to the upper band 
even though the gap vanishes. This is actually the same 
phenomenon as Klein tunneling for a ID version of the 
graphene bilayer, see e.g. the appendix of Ref. [22j . 
When in the gapless phase d < 0, this may be seen as 
two successive perfect Klein tunnelings for a ID mass- 
less Dirac electron: first going from the lower to the up- 
per band with unit probability and then going down to 
the lower band with certainty at the second Dirac cone. 
When k is non-zero but small, one can solve the coupled 
differential equations in perturbation theory as shown be- 
low, and show that the probability becomes finite. 
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B. Adiabatic basis 

It is also useful to write the same problem in the adi- 
abatic basis, which corresponds to diagonalizing H(t) 
with t being treated as a parameter. We call E a (t) = 
a \/ (t 2 + d) 2 + k 2 = aE + the adiabatic eigenenergies 
(plotted in Fig. [2^b) when d < 0, see also Fig. [I]), where 
a = ± is the band index, and \ip a (t)) the corresponding 
eigenvectors. They satisfy H(t)\ip a {t)) = E a (t)\ip a (t)} . 
The angle 6(t) is defined by s'md = k/E + and cos# = 
(t 2 + d) /E + , which allows us to write the adiabatic eigen- 
vectors as 

They form an orthonormal basis at each t. The state of 
the electron at any time can now be expressed in this 
basis as 

\m) = J2 A ^t)e- i - ftdt ' E ^ t ' ) \Mt)) (7) 

a 

in terms of two unknown amplitudes A a (t), which satisfy 
J2 a l^«l 2 = 1- The initial state is such that A_(— oo) — 
1 (up to a global phase factor) and A + (— oo) = and we 
are interested in P = |j4 + (oo)| 2 . Indeed, as t — > ±oo, 
9 w and \4>-(t)) w -|2) and |^+(*)) » Therefore 
at both initial and final times, the adiabatic and diabatic 
basis coincide. 

The time dependent amplitudes satisfy the following 
Schrodinger equations 

A + +A + (^ + ) = -{i> + \i>_)A_e i f tdt ' E +-W 
A_ + A_{il)_\i}_) = -(V^|^+)A +e -^ tti *' £ +-(*')(8) 

where E + - = E + - E-. As \ip±) = =f||^ t ), one has 

(ip±\ip±) = and = -{i/)+\4>-) = -f , so that 

the equations become 

A+ = -{i }+ \i ) .)A-e i ^ dtE +-^ 

A- = ((i/j + \iP_))*A+e- i f tdtE +-M (9) 

and are quite similar to the ones obtained in the diabatic 
basis, see eq. d5j). They also depend on two functions of 
time: one is the energy difference between the two basis 

states £_| (t) (instead of Ei 2 (t)) and the other is the 

coupling between these states (i() + \ip-)(t) =6/2 (instead 
o£iH 12 (t)). 

If the force is small, the motion of the electron is slow 
and k = c y p y (2m* ) 1/3 / (hF) 2/3 > 1 is large. This is 
the adiabatic limit and k can be thought of as an adi- 
abaticity parameter. It is also the semiclassical limit as 
it is equivalent to h — >• (in a purely classical problem, 
the transition probability would always be zero). In this 
limit, the electron stays in the state \tp-(t)), which in 
both limits t —> ±oo is ~ |2). As a consequence the 
transition probability P — |(l|^(co))| 2 « (oo)) | 2 



is also zero. When k is large but finite, it is possible to 
compute the transition probability in perturbation the- 
ory (this time the small parameter being 1/k) as shown 
below. 

To summarize, both when k 1 and /;> 1, the tran- 
sition probability vanishes. Away from these two limits, 
the probability will be non-zero. This already shows that 
the probability P is a non monotonic function of the per- 
pendicular gap k, which is in stark contrast to the LZ 
problem of a linear avoided band crossing. In the follow- 
ing, we use perturbation theory to compute the transition 
probability first in the diabatic and then in the adiabatic 
basis. 

III. STUCKELBERG THEORY IN THE 
GAPLESS PHASE 

We start by recalling the results we obtained previ- 
ously using the Stiickelberg theory in the gapless phase 
(d < 0), see the supplemental material of [T3]. We first 
compute the transition probability associated to the two 
successive LZ events, in the limit where they can be con- 
sidered to be far apart (i.e. deep in the gapless phase) 
using the Stiickelberg approach [T3HI5] . During a single 
LZ event the probability amplitude to stay in the up- 
per/lower band is ^/l — Pze TiVst where the Zener prob- 
ability is P z = e- 2lr5 . The non-adiabatic phase delay 
^ftpst (where =p refers to the upper/lower band) is given 
in terms of the Stokes phase [15] 

<p S t = ~ + 6(ln8-l) + Axgr(l-id) (10) 

where 

k 2 

5=^— (11) 

is the adiabaticity parameter in the linear band crossing 
problem (not to be confused with fc, which is the adia- 
baticity parameter in the quadratic band crossing). In 
the diabatic limit, 5 — > 0, the Stokes phase is ir/4 and it 
monotonically goes to zero in the adiabatic limit 5 — > oo. 
If the sequence between the two tunneling events is coher- 
ent, the two avoided linear crossings realize a Stiickelberg 
interferometer. The total probability amplitude to go 
from the lower to the upper band is the sum of the am- 
plitude for two distinct paths. In the first path, the elec- 
tron jumps to the upper band at the first Dirac cone and 
then stays in the upper band at the second, such that 
the amplitude is A + = -\/Pz x e iip + x - P z e~ lvst 
where —\[P~z is the amplitude to jump at the first avoided 

band crossing and tp + = J^~J—dtE + (t) is the phase dy- 

~ v 1^1 

namically acquired by the electron traveling in the up- 
per band between the two Dirac cones, with 2-y/jdf the 
time needed to travel between the two Dirac points. In 
the second path, the electron stays in the lower band 
at the first Dirac cone and then jumps to the upper 
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band at the second. The associated amplitude is A_ — 
yl - Pze l f st x e lip - x y/Pz where v^Pz is the ampli- 
tude to jump at the second avoided band crossing and 

(f- = Jv^L dtE— (t) is the dynamically acquired phase 
— v 1^1 

of the electron traveling in the lower band from one Dirac 
cone to the other. Note that the jumping amplitudes 
Ta/Pz at the two avoided crossings are opposite to each 
other. This is related to the fact that the linear LZ prob- 
lem is not exactly the same for the two Dirac cones: in- 
deed the local low-energy hamiltonians are slightly dif- 
ferent just as the ones describing the two different valleys 
of graphene [23]. The total probability P = \A + + A_\ 2 
is therefore 

P = 4P z (l-P z )sin 2 (^+<^ t ) (12) 

where (p,iy = tp— — ip+ is the dynamically accumulated 
phase between the two tunneling events [13] 

<Pdy = f^E+^dt = 4|d| 3 /2/( A) (is) 
J-yf\d\ \ d \ 

written in terms of the integral I(x) = 
j du\J (u 2 — l) 2 + x 2 . This probability as a func- 
tion of d and k is plotted in Fig. [3|a). 

If the two tunneling events are incoherent, the inter- 
ferences are washed out, sin 2 — > 1/2 and the probability 
becomes P incoh = 2P Z (1 - P z ) (see Fig. [§b)). 

In the limit 5 > 1, P z = exp(-27r<5) 0, 1 - P z 1, 
ip st ~ and ifd y ~ 4fcy|d[. Therefore 

rvk 2 

PR^eTvW S in 2 hk-sAdU (14) 

The probability goes to zero exponentially because of the 
large gap, as is usual for tunneling process. 

In the opposite limit 5 <C 1, P z -> 1, 1 — Pz « 2ttS -> 
0, Lp st ~ tt/4 and ip dy 8|d| 3 / 2 /3. Therefore 

P«2^4^sin 2 ^|d| 3 / 2 + ^ (15) 

VW\ V3 1 4; 

The probability also goes to zero but as k 2 because of 
the special symmetry when k — (conservation of the 
pseudo-spin a z ). 

Quantitatively, the Stuckelberg approach is valid if the 
Zener tunneling time ~ max((5, yS)/k is shorter than the 
time 2y|df it takes for an electron to travel between the 
two Dirac points |15j . This means that — d 1 and 
— d 3> k. Therefore, one needs to be deep in the gapless 
phase (far from the merging) and with a perpendicular 
gap that is not too large. 

IV. PERTURBATION THEORY IN THE 
DIABATIC/SUDDEN LIMIT 

In the diabatic basis, we perform perturbation the- 
ory in the perpendicular gap k <C 1. Assuming that 




FIG. 3: (Color online) Contour plot of the transition probabil- 
ity P computed with the Stuckelberg approach as a function 
of the merging gap d and the perpendicular gap k. The white 
region corresponds to d > 0, where the Stuckelberg approach 
is not defined, (a) In the coherent case, interferences both as a 
function of d and k are clearly visible as well as the vanishing 
of P in the k — > and k — > 00 limits. The probability varies 
between and 1 as given by the color code (color steps corre- 
sponds to 0.1). (b) In the incoherent case, the oscillations are 
washed out and the maximum probability is 1/2 instead of 1 
in the coherent case (note the change of color scale for P). 



A 2 (t) rj 1 for all t gives the probability P = |Ai(oo)| 2 to 
tunnel from the lower to the upper band in terms of the 
amplitude 



Ai(oo) = -k I dtA 2 (t)exp[i / dt'E 12 {t')} 

J ~oc JO 

/CO 
dtexp[i(2td + 2t 3 /3)} (16) 
-OO 
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computed at first order in fc. The probability 

P« 4 2 / 3 7r 2 fc 2 [Ai(4 1 / 3 d)] 2 (*<1) (17) 

is given in terms of the Airy function which has the fol- 
lowing definition (when its argument x is real): 



Ai(ar) 



1 

2^ 



dye 



i(ly 3 +xy) 



(18) 



As it will be useful later, we also perform a saddle point 



analysis of the integral in eq. ( 16 ) in three different limits 
to obtain simpler analytical results. If d =^ 0, there are 
two saddle points t in the complex time plane. If d > 0, 
to = iiVd and only to — i\fd contributes, as Im to > 
is needed. If d < 0, to = ±\/Pl anc ^ the two saddle 
points contribute giving rise to interferences (this is really 
a stationary phase approximation). If d = 0, there is a 
single saddle point at i = 0. The results of the saddle 
point approximation are 




P 



P 



4/3 , : 



fc 2 if \d\ < 1 (19) 



FIG. 4: (Color online) Contour plot of the transition probabil- 
ity P as a function of the merging gap d and the perpendicular 
gap k. For small k, it is computed with the diabatic pertur- 
bation theory eq. (171; while for lar ge fc , it is computed using 
adiabatic perturbation theory eq. (30 1, see section V. Inter- 



P 



2nk 2 



TTfc 2 

2^fd 



Sd 3 / 2 



if d> 1 



3/2 



if - d > 1 



(20) 



(21) 



ferences as a function of d are clearly visible in the gapless 
phase (d < 0). In the gapped phase (d > 0), the probabil- 
ity vanishes exponentially. White regions corresponds to the 
probability exceeding 1. Indeed close to k ~ 1, the two per- 
turbative approaches break down. The color code is the same 
as in Fig. [3ja) 



The last case recovers the result of the previous section, 

see Eq. featuring Stiickelberg oscillations. These coupling 9 is small, we find that A_(t) « 1, A 



with the initial conditions A_(— oo) = 1 (up to a phase 
factor) and A + (— oo) = 0. If we now assume that the 



three limits are well known expansions of the Airy func- 
tion. The transition probability P of eq. (171 is plotted 
in Figs. [4] and [5] It goes to zero as k — > as expected 
from the a z conservation and increases quadratically with 
k. As a function of d, P shows oscillations when d < 0, 
which we interpret as Stiickelberg interferences, and de- 
creases exponentially when there is a true gap d > in 
the diabatic spectrum. 



V. PERTURBATION THEORY IN THE 
ADIABATIC /SEMICLASSIC AL LIMIT 

In order to consider the opposite adiabatic limit (fc 
1), we use perturbation theory in the adiabatic basis |24j . 
From the adiabatic eigenenergies and eigenvectors, we 
find that £+_ = 2E + , = and {ip+\ip-) = 

8/2 = —kt/E^. The Schrodinger equations to be solved 
arc therefore 



.| e 2i/*dtB+(t) and therefore: 



A+(co) a - 



-|-OG 



dt-e 



2i p dt'E+(t') — 



dt-e 



i<t>(t) 



A_ = -A + e-~ 2i ^ dt ' E + ( - t " > 



(22) 
(23) 



(24) 

where 4>(t) = 2 J dt'E + (t') is the adiabatic phase and 
ti is the lower bound of the phase integral, which is un- 
decided for the moment except that it has to be a real 
number. 

The amplitude ^4+(oo) can be computed by integra- 
tion in the complex plane. Firstly, there are four poles 
corresponding to E + (t) = i.e. to t 2 = —d± ik. Note 
that these are also saddle points as cj>(t) — 2E + . In the 
following we refer to them simply as poles even when they 
are playing there the role of saddle points. If we write 
— d + ik = \fk 2 + d 2 e l P , which defines the angle /?, the 
four poles are h = (fc 2 + d 2 f^e lfi l 2 , t 2 = t\, t 3 = -t x 
and £4 = — t*. In addition to the four poles, there are 
also branch cuts coming from the square root function 
in the exponential. The corresponding branching points 
are at the same position as the poles. Therefore, there 
is a branch cut linking t\ and t<± and another one link- 
ing ti to £3, for example. When constructing a closed 
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FIG. 5: Transition probability P as a function of d for fixed 
k = 0.1 as computed in diabatic perturbation theory, see 
cq. (17|. For negative argument (gapless phase, d < 0) it 



shows Stiickelberg oscillations and then decays exponentially 
for positive argument (gapped phase, d > 0). There is an 
inflexion point right at the merging (when the argument van- 
ishes d — 0). 



contour, one has to keep in mind this branch cut struc- 
ture. Of the four poles, only t\ and £4 have a positive 
imaginary part and are therefore relevant as we want to 
close the integration contour in the upper plane (see Fig. 
[6]). Since Im t\ = Im £4, both poles contribute equally to 
the amplitude. We choose the lower bound ti = such 
that 4>{t) — 2 L dt'E + (t'). It is important to make a sin- 
gle choice for t; for both poles as they will interfere. As 
E + (-t) = E + (t) and t A = -t\, we have cj){U) = -cj)(ti)*. 
One possibility is therefore to construct a contour that 
encloses both these two poles and the branch cuts. 




FIG. 6: (Color online) Poles ti and £4 = — ti in the complex 
time plane. Poles are represented at the merging transition 
(d — such that arg t 1 = 7r/4 = — arg £4). The short green 
arrows indicate their motion when d increases (gapped phase 
d > 0) and the long red arrows that when it decreases (gapless 
phase d < 0). When the poles have a finite real part, there 
are oscillations in the probability; whereas a finite imaginary 
part implies an exponential decay of the probability. Note 
that there is a branch cut relating the two poles. 

Another trick to perform this integral is to make a 
change of variable from the time t to the phase </> variable, 
see for example Ref. [27]. The resulting integral 



1 r + °° 



dO 

6 

d<j> 



,i<f> 



(25) 



is over a function [d6 / 'd<\>)e % ^ that has no branch cut 
anymore and only four isolated poles at <pi = 4>{t\) 1 
4>(t 2 ) = 0i, (f>(t 3 ) = -4>i and cj)(t 4 ) = —(f>*. The residue 
theorem can now be used with a simple contour closed 
in the upper complex plane of the variable. As the 
residues are — e 1 ^ 1 /3i and e -4 ^ 1 /3i, we obtain 



A+(oo)«^(e^-e- 



2l7T 



sin(Re <fii)e~ 



-Im < 



This is valid whatever the sign of d. Therefore 



4t: 2 



P ps — sin 2 (Re <t>i)^ 2lvci 01 



(26) 



(27) 



where 



0i = (p{ti) = 2k z ' 2 / du^l + iu 2 + D) 2 (28) 
Jo 

with D = d/k and ui = ti/Vk = {\f^/l + D 2 - D + 

iVVl + D 2 ~+D)/V2. 

The result we found is the first-order perturbation in 
the adiabatic basis. The exponential behavior is cor- 
rect but not the prefactor, even in the adiabatic limit, 
as argued by Landau long ago [2SJ [25] . This is known 
in the literature as the "7r/3 problem" [27J [2S]. Actu- 
ally in the case of a single linear band crossing {E\ = 
—E2 = at/2 and H12 = constant), which is the stan- 
dard LZ problem, adiabatic perturbation theory gives 
P = (7r/3) 2 exp(-27r|i?i 2 | 2 /a) US whereas the exact re- 
sult found by Zener is P z = exp(-2n\H 12 \ 2 /a) [T7J. The 
reason for this discrepancy is well explained in [30) : it 
comes from the fact that each order of the adiabatic per- 
turbation expansion for A + contains a term of the form 
#exp(— n\Hi2\ 2 /a). Obtaining the exact factor in front 
of the exponential requires re-summing the whole series 
by keeping only the dominant - in the adiabatic limit 
exponential behavior in each order. This series has a first 
term which is it /3 and a sum which is 1 [24] [28H30] . In 
the adiabatic limit, the correct pre-exponential factor in 
the usual LZ problem is such that 



P : 



e = e 



21 



(29) 



which amounts to reducing the residue of the pole found 
in first order adiabatic perturbation theory from 7r/3 to 1 
in the amplitude. This can also be done in the two poles 



case and we find that A + (po) 
e i<t>{t\) 

of eq. 



(tl) 



i<K*i)* 



(27), 



l ' so that the probability becomes, instead 



P«4sin 2 (Re <j>i)e 



-2Im ■ 



(28) 



(k > 1) (30) 
As we will later see, 



where (f>i{k, d) is given in eq. 
this result agrees very well with the exact numerical solu- 
tion. It also agrees with the Stiickelberg theory in the adi- 
abatic limit when Pz <C 1 such that P » APz sin 2 (y>dy/2) 
where Pz = cxp(— ixk 2 j{2\J— d)) deep in the gapless 



phase (indeed 2e 



-21m . 



). Therefore, we 



take this result as the correct analytical expression in the 
adiabatic limit. 

We now come back to the phase <f>i given in eq. [28] 
As ui{D) = iui(-D)*, one has 0i(fc, d) = i<j>i(k,-d)* 
and therefore Re 0i(fc, d) = Im </>i(fc,— d), which allows 
one to express P in terms of Im^i only. The integral 
J(D) = Im /J* 1 duy/l + (u 2 + Df giving Im 0i(fc,d) = 
2fc 3 / 2 J(d/fc) can be computed numerically for any D and 
analytically in three limits. When D ~ 0, J(D) « 



r(i/4) 2 



+ 



r£>. When D ->■ oo, J(D) « 2L> 3 / 2 /3 



12^ ^ T(l/4) 2 

]n D / (4y/~D) . In practise, a good approximate interpola- 
tion between and oo is J(D) w + §L> 3 / 2 . When 

D -> -oo, J(D) w 7r/(8v^D|). This function J(D) is 
plotted in figure [7] 




FIG. 



7: 



(Color online) Integral J = 
Im f Q ul du^/l + (u 2 + D) 2 - giving Im 0i = 2k 3/2 J(d/k) 
- plotted as a function of D. The numerical calculation is 
in continuous red and is compared to different analytical 

r(i/4) 2 , 



tD in 



results: 7r/(8y|D|) in dashed green, [ 2 r~ — ^ 



1 T(l/4)2 

dotted magenta, and 2D 3/2 /3 + ln£)/(4V f D) in dot-dashed 
blue. Interpolation formula (for positive d) — h |/) 3//2 
is in thin black. 

From the behavior of J(D), we can obtain approximate 
analytical results for the probability P in three limits: 



P « 4e-^ fc3/2 sin 2 (^M!fc 3 / 2 ) if |d, « k (31 ) 

6Vf 



P«4 e -^sin 2 (^l) if d> fc > (32) 
4Vc? 



Psa4e 2 \^sin 2 ( 
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FIG. 8: (Color online) Transition probability P as a func- 
tion of d (at fixed k = 2.5) computed in adiabatic perturba- 
tion theory (continuous red line), see eq. (30 1. Also shown, 
in dashed black, is the incoherent probability Pi nco h, see eq. 
(341. 



if - d > k > 



(33) 



It is interesting to discuss the motion of the two poles 

h = (k 2 + d 2 ) 1 /4 e «/3/2 and ti = _^ where _ d + ik = 

Vk 2 + d 2 e 1 ^ , in the complex time plane as k and d vary 
[3"T] . These two poles correspond to the band crossings 
at complex times and always exist (whatever the sign of 
d and even at the merging). When d — 0, [5/2 = 7r/4 
and Re t x = Im t x . When d > 0, tt/4 < [5/2 < vr/2 
and Re t\ < Im t±, the poles are close to the imaginary 
axis, the corresponding exponentials are essentially de- 
caying and the probability as well. In the limit d — > +oo, 
[5/2 — > n/2 and the two poles are on the imaginary 
axis. Remember that, in the diabatic limit k — ¥ 0, we 
found a saddle point at to = i\fd, i.e. (5/2 ~ 7r/2 and 
(fc 2 + d 2 ) 1 / 4 ~ Vd. When d < 0, < [5/2 < tt/4 and 
Re t\ > Im t\, the poles are close to the real axis, the 
corresponding exponentials are essentially oscillating and 
the interference of the two give oscillations in the prob- 
ability. In the limit d —> — oo, [5/2 — > and the two 
poles are on the real axis. Remember that, in the di- 
abatic limit k 0, we found two stationary points at 
t Q = ±V^d, i-e. ,8/2 ~ and (fc 2 + d 2 ) 1 / 4 - V^d. The 
approximate Stiickelberg theory also falls in this general 
frame. It corresponds to a situation where — d 3> 1, fc 
such that t\ ps \/\d\ and ti ~ — yfdf (there, we identified 
ti — ti w 2-v/fd[ as the time needed to travel between the 
two Dirac cones). We speculate that in the general case 
of arbitrary k and d, there are always two separated com- 
plex poles with the same positive imaginary part. The 
motion of the poles in the complex t plane as d changes 
at fixed fc ^ is illustrated in Fig. |6j 



If the contribution of the two poles add incoherently, 
sin 2 — > 1/2, and the oscillations are washed out: 

Pncoh « 2e- 2Im ^ (34) 

as in the Stiickelberg theory when Pz <C 1. 



VI. MODIFIED STUCKELBERG FORMULA 

In the preceding section, adiabatic perturbation the- 
ory helped us uncover a general two poles structure - 
either in the complex t or complex plane -, which 
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leads to a total probability of the Stiickelberg form 
P = 4Ps ( 1 — Ps ) sin 2 ( . . . ) , where Ps is the probability for 
a single avoided crossing. This should be valid for all k 
and d and not only when the spectrum is gapless. A rea- 
sonable guess (see also ;32J) is to combine the adiabatic 
perturbation theory, giving the exponential weight of the 
two poles and their interferences in the adiabatic limit, 
with the Stiickelberg approach, giving the Ps(l — Ps) 
structure. From eqs. (12) and (30), we obtain: 



P«4e 



-21m . 



*(1 



-21m . 



L )sin 2 (Re0i + ^„ Q ) (35) 



There, e 4< " is the amplitude to tunnel for a single pole, 
so that e~ 2Im ^ 1 plays the role of the Zener probabil- 
ity Pz for a single Dirac cone and Re <j>\ that of (pd y /2. 
The quantity (p na is the non-adiabatic phase acquired 
by a particle when it does not tunnel at a single pole 
- the associated amplitude being y/l — e _2Im 0ie J¥, ™ a . 
We only know its expression in the Stiickelberg limit 
(d <C — 1, —k), where it is given by the Stokes phase tpsti 



see eq. (10) with 8 as in eq. (11). Here, we assume that 
f na ~ VSt f° r au ^ an d d, which is a reasonable approx- 



imation except when k < 1 and d > 0. Equation (35) 



should be exact both for small k and negative d, where 
it recovers the Stiickelberg result eq. (12), and for large 
k, where it recovers the result of adiabatic perturbation 
theory eq. (30) for all d. By continuity, it should also 
be reasonable in the intermediate region k ~ 1, see Fig. 
13 b). It allows one to have an approximate analytical 



formula that can describe the crossover from small to 
large k at fixed d, see Fig. [12] This modified Stiickelberg 
probability is plotted in Fig. [9]ja). As seen, this for- 
mula is not applicable for positive d and small k as the 
relevant non-adiabatic phase is no more simply given by 
the Stokes phase. In the incoherent case, the probability 
becomes 




P, 



2e 



-21m 



(1 



-21m 



11 ) 



(36) 



and is plotted in Fig. ^b). As this incoherent probability 
does not depend on the partly unknown non-adiabatic 
phase, it should be reasonable in the whole (d, k) plane. 



VII. NUMERICAL SOLUTION AND 
COMPARISON BETWEEN DIFFERENT 
APPROACHES 



FIG. 9: (Color online) Contour plot of the modified 
Stiickelberg transition probability P as a function of the merg- 
ing gap d and the perpendicular gap k. (a) Coherent case (see 
eq. (35)): the probability is in between and 1 (the color 
code is the same as in Fig. [3^a)). Note that the modified 
Stiickelberg formula does not work in the (d > 0, k < 1) re- 
gion as the non-adiabatic phase is not properly given by the 
Stokes phase; (b) Incoherent case (see eq. ( |36[ )): the proba- 
bility is in between and 0.5 (the color code is the same as 
in Fig. §b)). 



The coupled first-order differential equations of sec- 
tion II, see eq. ^ and eq. ([9]), are solved numerically. 
We checked that solving these equations either in the 
diabatic or in the adiabatic formulation gives the same 
answer (up to numerical errors of order 10 -3 in the prob- 
ability). We can therefore consider that these numeri- 
cal solutions are essentially exact and use them to check 
the approximate analytical solutions. The probability 
obtained numerically for any d and k is shown in Fig. 
[TTJ] When compared to diabatic perturbation theory, the 
agreement is perfect for small k <C 1. When compared 



with the Stiickelberg theory, the agreement is very good 
when d is very negative and k not too large compared to 
—d (—d ^> 1 and —d^k). It also compares very well 
with adiabatic perturbation theory (provided 7r/3 — > 1) 
when k is large (k ^> 1). 

To compare the different approaches, we first concen- 
trate on the d — case exactly at the merging transition. 
The numerical solution along with the diabatic and adi- 
abatic perturbative results are shown in Fig. [TT] Note 
the excellent agreement in both the k — > (diabatic per- 
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FIG. 10: (Color online) Contour plot of the numerically com- 
puted transition probability P as a function of the merging 
gap d and the perpendicular gap k. Oscillations are clearly 
visible in the gapless phase, whereas the probability is van- 
ishingly small in the gapped phase. The vanishing of the 
probability in both the diabatic < 1 and adiabatic k 3> 1 
limits is also visible. The color code is the same as in Fig. 
§a). 



turbation theory) and the k — > oo limits (adiabatic per- 
turbation theory). At large k, a surprising oscillation 
in the probability is seen both in the numerical solution 
and in the adiabatic perturbative result. It is surpris- 
ing because the spectrum (whether diabatic E = ±t 2 or 
adiabatic E = ±^Jt 2 + k 2 ) features at most a single real 
time crossing. However in complex time, the adiabatic 
bands cross twice. This leads to an interference between 
the two complex poles t\ and t± and results in oscillations 
in the probability P. 

Next we consider the gapless region d = — 1 and com- 
pare the different analytical approaches to numerics, see 
Fig. [12] Note the excellent job done by the modified 
Stiickelberg formula which is able to describe the whole 
crossover from small to large k. The only discrepancy 
with the numerical result is close to the maximum prob- 
ability near k = 0.5. 

Then we consider small k and compare numerics, the 
Stiickelberg approach and diabatic perturbation theory 
as a function of d, see Fig. 



theory agrees very well wit 



13 a) . Diabatic perturbation 



r the numerical result except 



for a small difference close to d = — 1 where the probabil- 
ity is not small and the approximation is therefore not so 
good anymore. Stiickelberg theory works very well deep 
in the gapless phase and its validity breaks down as one 
approaches the merging transition. The opposite limit 
of large k shows that the adiabatic perturbation theory 
is very good (see Fig. [l3[c)). The Stiickelberg theory 
works qualitatively in the gapless regime but not as well 
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FIG. 11: (Color online) Transition probability P at the merg 
ing d = as a function of k. The numerically exact result is 
in continuous blue. The diabatic perturbative result eq. ( 17 
is in dashed red and the adiabatic perturbative result eq. ( 30 



is in dotted black, (a) k between and 2. (b) k between 1.5 
and 3 (note the change of vertical scale by a factor 10 3 ): there 
is a tiny oscillation due to the interference between two poles. 



as for small k. There are also regimes which are difficult 
to access analytically. This is the case for intermediate 
k ru 1. See Fig. 13 'b) for the k = 1 curve as a function of 
d. Stiickelberg theory works fine but only covers the large 



— d regime, whereas both perturbative calculations (not 
shown for clarity) are unreliable for intermediate k. The 
modified Stiickelberg is qualitatively correct also when 
fc~ 1. 



VIII. COMPARISON TO THE EXPERIMENT: 
ABSENCE OF INTERFERENCES 

Recently, an experiment with ultracold fcrmionic 
atoms in an optical lattice could study the merging tran- 
sition and the inter-band tunneling of atoms performing 
Bloch oscillations [BJ. There, atoms moving in an arti- 
ficial graphene-like crystal could mimic Bloch electrons 
in a usual solid state crystal. In Ref. [13], in order to 
understand the result of this experiment, we computed 
the inter-band transition probability for a single atom 
using the approximate Stiickelberg theory as a function 
of k and d. Then we translated these two parameters 
into the experimentally tunable laser intensities V x and 
Vx defined in [BJ . Qualitatively, V x controls the merging 
transition and is roughly equivalent to — d (called —A* 
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(a) 
k=0.1 



FIG. 12: (Color online) Transition probability P at fixed d = 
— 1 as a function of k. The numerically exact result is in 
continuous blue, the diabatic perturbative resu lt eq. \17\ is 
in dashed red, the Stiickelberg result eq. ( 12 I is in dotted 
black, the adiabatic perturbative result eq. ( 30 1 is in dot- 



dashed green and the modified Stiickelberg formula eq. 
is in long dashed magenta. 



in [T3]), whereas Vx controls the transverse gap and is 
equivalent to k (called c x in [13]). A last step was to 
average the probability over the atomic distribution of a 
two-dimensional degenerate Fermi gas. The agreement 
between theory and experiment was found to be very 
good: compare Fig. 4(b) in [6] with Fig. 4(b) in [13] . 

However, as the Stiickelberg theory is only valid in the 
gapless phase (d < 0) and not too close to the merging 
transition (see Fig. 3[b)), we could only compare theory 
and experiment in the gapless region. Near the transition 
and in the gapped region, the experimental signal was 
vanishingly small and could not be compared with any 
theoretical prediction. Within the present framework, it 
is now possible to understand the inter-band probability 
very close to the merging. When looking at Fig. 4(b) 
of Ref. 6J in detail, one sees that the red line of maxi- 
mum transition probability - which lies essentially in the 
gapless region - actually crosses the merging line and 
slightly extends in the gapped region at very small Vx- 
This line qualitatively corresponds to Pz = 1/2 such that 
P = 1/2. Such a behavior is found in our calculations as 
well: see Fig. [9^b), where the orange region of maximum 
probability (between 0.4 and 0.5) lies essentially in the 
gapless (d < 0) region but also slightly extends to the 
gapped region (d > 0) reaching d ~ 0.4 when k — > 0. 

We now consider the merging point (d = 0) and study 
both the inter-band probability P 2 for the motion in the 
direction where two Dirac cones are hit (x direction) and 
that, Pi, in the perpendicular direction in which a single 
cone is hit (y direction) [19) . In the present paper, we 
concentrate on P 2 - which is called P(d, k) - as Pi is 
simply given by the usual LZ formula and was studied 
in detail in Ref. [13) . The merging point is special in 
the sense that the spectrum is gapless and features a sin- 
gle contact point, which is linearly dispersing in the p y 
direction and quadratically in the p x direction [S]. The 




FIG. 13: (Color online) Transition probability P at fixed k as 
a function of d. The numerically exact result is in continuous 
blue, the diabatic perturb ative result eq. ( |17| in dashed red, 
the Stiickelberg result eq. ( 12 1 in dotted black, the adiabatic 
perturbative result eq. ( |30[ ) in dot-dashed green and the mod- 
ified Stiickelberg formula eq. ( |35[ ) in long dashed magenta, 
(a) k = 0.1 (the modified Stiickelberg formula coincides with 
the Stiickelberg probability when d < - it is not shown for 
clarity - and is not applicable when d > 0); (b) k = 1; (c) 
k = 2.5. Note the different probability scales in the three 
graphs. 



LZ formula gives Pi = exp[— -K{p 2 / (2m*)) 2 / (hc y F)]. In 
the coherent case, the numerical solution of section VII 
gives P 2 = P(d = 0,k) where k = c y p y (2m*) 1 / 3 /(hF) 2 / 3 
(see the continuous blue line) and, 
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as plotted in Fig. 
in the incoherent case, the modified Stiickelberg eq. ( |36[ ) 

gives P 2 



-^^[l - exp(— 

The probability Pi depends on p x and varies between 



2exp( 



12 



and 1, whereas P2 depends on p y (i.e. on k) and varies 
between and ~ 0.55 (coherent case, see Fig. [TT|a)) 
or 0.5 (incoherent case). The ratio P2/P1 can therefore 
take any positive value depending on what are the rel- 
evant p x and p y values. The latter depend on the size 
of the atomic cloud and on the way the averaging over 
the atomic cloud is done. For example, for a single atom 
Px — Py = giving Px — 1 and P 2 = so that P2/P1 = 0. 
In particular, there is no reason for this ratio to take the 
simple value 0.5 [6]. We have performed averaging over 
various atomic cloud sizes comparable to that in the ETH 
Zurich experiment and find that (P2} / (Pi) can vary be- 
tween and ~ 0.7. 

One very striking experimental fact remains to be ex- 
plained: the agreement is actually obtained with the 
incoherent inter-band probability (see e.g. Fig. [9](b) ) 
rather than with the coherent probability (see e.g. Fig. 
10 1. In other words, Stiickelberg oscillations (interfer- 
ences) are not observed in the experiment, whereas they 
are predicted. Here we would like to discuss this specific 
point in more details. The absence of interferences could 
be due to (i) decoherence, (ii) blurring or (hi) washing 
out because of some averaging process, (i) Decoherence 
is unlikely in a cold atom experiment with almost non- 
interacting fermions. We estimate the decoherence time 
due to spontaneous emission following Ref. [33 . It is 
roughly given by 1/7 = (<5/Sl) 2 /7 ~ 10 3 s where 7^6 
MHz is the natural line width for the relevant transition 
of 40 K, 8 ~ 108 THz is the detuning and - 1 GHz is 
the Rabi frequency estimated from M7 2 /S ~ Er, where 
En r~j 4.4 kHz is the recoil energy. It is much longer than 
the experimental time, therefore ruling out decoherence 
as a possible mechanism to explain the absence of in- 
terferences, (ii) Blurring of the interferences could also 
occur because of the detection process using a finite pixel 
size. We checked that possibility and found that the pixel 
size is small enough that it should allow experimentalists 
to resolve the interferences, (hi) We are left with the 
possibility of washing out of the oscillations due to sev- 
eral averaging processes. We included averaging over a 
two-dimensional atomic distribution in reciprocal space, 
which only resulted in slightly smoothing the oscillations 
(compare Fig. 10 here and Fig. 4(d) in [H]). How- 
ever, the atomic cloud in the experiment was actually 
not two but three-dimensional, even though the optical 
lattice was two-dimensional. The atomic gas was indeed 
confined by an anisotropic three dimensional harmonic 
trap but very far from the regime where one of the di- 
rection of motion would be frozen. This means that the 
system is best seen as a bunch of parallel one-dimensional 
tubes, each tube corresponding to a single site of a two- 
dimensional honeycomb-like lattice. The atoms hop in a 
kind of tight-binding lattice in the xy plane (except for 
a weak harmonic trap muj 2 x 2 /2 + muj 2 y 2 /2) and are al- 
most free to move in the z direction (except for a weak 
harmonic trap muj 2 z 2 /2). The period of the harmonic 
motion in the z direction 2tt/uj z is very long compared to 
the time an atom spends in the Stiickelberg interferome- 



ter ~ 2y/\d\tca r , where t car — (2m*ft) 1 / 3 /F 2 / 3 . One can 
therefore think that an atom moves in the interferometer 
at an almost constant z. However, because of the finite 
laser waist, the laser intensities are inhomogeneous, so 
that atoms at different z experience a slightly different 
optical lattice. In other words, the parameters d and k 
of the universal hamiltonian are slightly z-dependent. As 
seen in Fig. [TUJfor example, the interferences in the inter- 
band probability P are essentially a function of d (and 
not so much of k) , with a fringe spacing of roughly 5d ~ 1 
(which is the same as saying that <JA* ~ O.OAEr |13j). 
From the experimental conditions of Ref. [5] , we estimate 
a laser waist of ~ 150 microns and a cloud radius of ~ 30 
microns in the z direction (half of the tube's length) so 
that the parameter d varies by roughly 0.7 between the 
center and the edge of the atomic cloud. As this is com- 
parable to the spacing between a dark and a bright fringe, 
it should be enough to wash out the oscillations. In the 
experiment, the inter-band probability is automatically 
averaged over the third spatial direction, i.e. along the 
tubes axis. Therefore, we think that the averaging over 
the third spatial direction could be responsible for the 
absence of the oscillations in the inter-band probability. 
An alternative explanation for the disappearance of the 
oscillations was very recently proposed in Ref. [34] . It is 
based on the spatial inhomogcncity of the applied force in 
the 2D plane, which also leads to averaging and washing 
out of the probability fringes. 

By breaking the inversion symmetry of the lattice, it 
is also possible to induce a mass to the Dirac fermions, 
i.e. to gap the Dirac cones when d < 6J. Such a sit- 
uation is easily incorporated in our theory by a simple 
mapping k — ¥ \Jk 2 + g 2 . The hamilonian (I3p becomes 
H = [t 2 + d\a z + ka x + go~ y where g is the (dimension- 
less) mass gap induced by inversion symmetry break- 
ing. The inter-band transition probability V(d, k,g) 
when g ^ is simply related to that P(d, k) at g = 
by V(d,k,g) = V(d,^WT^,0) = P(d,ypTF). 
This mapping is easily found by looking at the coupled 
differential equat ions ([5] ), in which H12 = k becomes 
H12 = k — ig~ y/k 2 + g 2 e tl where 7 = Arg (k—ig). The 
phase 7 is time independent and can be gauged away so 
that only the modulus of k — ig matters and _ff 12 becomes 

V k2 +9 2 - 



IX. CONCLUSION 

Inspired by a recent experiment probing the merging 
transition of Dirac cones via Bloch-Zener oscillations of 
ultracold fermionic atoms 1 1 3J , we have studied inter- 
band tunneling for a quadratic band crossing. The lat- 
ter problem depends on two dimensionless parameters, 
which are the merging gap d and the perpendicular gap 
k. We computed the probability P for a particle to tun- 
nel from the lower to the upper band as a function of 
k and d. Qualitatively, the probability oscillates as a 
function of d in the gapless phase and decays exponen- 



13 



tially in the gapped phase. The oscillations are a result of 
Stiickelberg interferences. As a function of k, the prob- 
ability shows quite an unusual non-monotonic behavior: 
P vanishes exponentially in the adiabatic/semiclassical 
limit (large k) , which is the expected tunneling behavior 
in the large gap limit, but it vanishes also in the opposite 
diabatic/sudden limit (small k) as a result of a special 
symmetry. Indeed, when k = 0, the conservation of the 
pseudo-spin a z implies that P vanishes. When k ^ 0, 
this symmetry is broken and, quite counter-intuitively, 
opening of a gap leads first to a quadratic increase of the 
probability to tunnel between the bands. In addition, 
when k ^> 1, there are oscillations of P (as a function 
of both k and d) whatever the sign of d. These are due 
to interferences between two poles in the complex time 
plane. The latter exist not only in the presence of Dirac 
points (gapless phase) but also in the gapped phase (in 
which case the bands do cross but at times with a finite 
imaginary part). 

The probability P of inter-band tunneling was calcu- 
lated using different methods. To summarize: the numer- 
ically exact solution of the time-dependent Schrodinger 
equation is given in section VI. We also used approximate 
analytical techniques to compute P: for small k <C 1 and 
arbitrary d, we used diabatic perturbation theory, see eq. 



(17 1. For negative d <C — 1 and small k <C — d, we em- 
ployed the Stiickelberg approach, see eq. (12). And for 



large k ^> 1, we used adiabatic perturbation theory, see 
eq. (30). For intermediate fc's, we have no exact analyt- 



ical prediction but an approximate modified Stiickelberg 
formula, see eq. (35), that compares well to the numerics 



in the whole negative d region and also for large k and 
positive d (adiabatic regime). Using the tools we have de- 
veloped, it should be possible to compute the inter-band 
tunneling probability for many two-bands hamiltonians. 

Note added: After completion of the present work, we 
became aware of related articles in the context of atomic 
collisions, in which a parabolic level crossing problem was 
studied, see Ref . [32] • The specific case d = (exactly at 
the merging transition) was also very recently analyzed 
in [35], where it is called parabolic level glancing. 
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